Noise at the Wigner Glass Transition 
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Using a simple model for interacting electrons in two dimensions with random disorder, we show 
that a crossover from a Wigner liquid to a Wigner glass occurs as a function of charge density. The 
noise power increases strongly at the transition and the characteristics of the 1/f" noise change. 
When the temperature is increased, the noise power decreases. We compare these results with recent 
noise measurements in systems with two-dimensional metal-insulator transitions. 
PACS numbers: 71.30.-Hh,73.40.Qv 



In the two-dimensional (2D) metal-insulator transition 
(MIT) regime, both the Coulomb interactions between 
electrons and the disorder are expected to be strong, lead- 
ing to the formation of an electron glass [1,2]. Recent ex- 
periments in 2D electron systems have revealed changes 
in the characteristics and the amplitude of the conduction 
noise as the change density of the system is varied from 
a high charge density metallic phase to a lower charge 
density insulating phase [3-5] . This has been interpreted 
as evidence for an onset of glassy dynamics near the in- 
sulating phase. These studies also find that the noise has 
a 1//" characteristic, with a — 1.0 in the metallic phase 
changing over to a ~ 1.8 in the glassy phase. Other ex- 
periments have found similar results with a ~ 0.75 in 
the metallic phase and a — 1.3 near the insulating phase 
[6]. The glassy phase is signified by a large increase in 
the noise power along with the change in a [3,4]. In 
addition, the noise power was observed to decrease with 
temperature [4,5], in contrast to single electron models 
with thermally activated trapping [7] and other models 
[8] that predict an increase in the noise power with T. 
This suggests the importance of electron-electron inter- 
actions at the MIT. Other recent experiments near the 
2D MIT have also found 1//" noise, strong increases in 
noise power with decreasing charge density, and decreas- 
ing noise with increasing T [6,9]. Additionally, 1//" noise 
fluctuations in thin granular films have been interpreted 
as evidence for a glassy electron state [10]. 

In theoretical studies, it was proposed that at the 2D 
MIT a freezing from an electron liquid to a partially or- 
dered Wigner glass [11] or a more strongly disordered 
electron glass [12,13] may occur. Other theories sug- 
gest that an intermediate metallic glass phase appears 
between the liquid and insulating phases [14]. It is also 
possible that the metallic glass phase may consist of solid 
phase insulating regions coexisting with string-like liquid 
regions. Studies of glassy systems often find coopera- 
tive string-like motions or dynamical heterogeneities [15]. 
Such motions can give rise to correlated dynamics and 
large fluctuations near a glass transition. It is, however, 
unclear what the origin of such cooperativity would be 
in the electron glass systems. 

The strong enhancement of the noise observed when 
passing from a liquid to a glassy phase in the 2D charge 



system has noticeable similarities to the noise found in 
studies of vortex matter in superconductors [16]. In the 
vortex system, the existence of a transition from a weakly 
pinned liquid phase to a more strongly pinned glassy 
phase has been established. Noise studies in low tem- 
perature superconductors have shown that, as the glassy 
phase is approached from the liquid state, the voltage 
noise has a 1//" characteristic with increasing a, corre- 
sponding to an increase in the noise power [16]. 

The vortex system studies suggest that similar physics 
may be occurring in the 2D electron systems when there 
is a transition from a liquid like phase with low noise 
power to a pinned glassy phase with high noise power. 
This is also consistent with the theoretical prediction that 
the electron liquid freezes into a 2D disordered solid. In 
this work we propose a simple model for a classical 2D 
electron system consisting of interacting electrons with 
random disorder and temperature. We monitor the fluc- 
tuations and noise characteristics of the current as a func- 
tion of electron density or temperature. The advantage 
of our model is that a large number of interacting elec- 
trons can be conveniently simulated, while a full quantum 
mechanical model of similar size would be computation- 
ally prohibitive. Despite the limitations of this model, we 
show that this approach captures many of the key experi- 
mental observations. Additionally, although our primary 
focus is to gain insight into the physics near the 2D MIT, 
our model is also relevant for other classical charge sys- 
tems undergoing transitions from glass to liquid states, 
such as charged colloids interacting with random disor- 
der. Our model is similar to previous studies of 2D classi- 
cal electron systems with disorder [17,18]; however, these 
previous studies focused on the microscopies of the de- 
fects in the lattice [17] or the sliding dynamics [18]. In 
the present work we focus on the noise fluctuations in the 
strongly disordered phase as it changes from a liquid to 
a frozen state as a function of electron density for a flxed 
amount of disorder. 

Our model consists of a 2D system of Ns interacting 
electrons with periodic boundary conditions in the x and 
y-dircctions. There are also Np defect sites which attract 
the electrons. We assume the electron motion is at fi- 
nite temperature and the time evolution occurs through 
Langevin dynamics. The damping on the electrons comes 
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FIG. 1. The average electron velocity v vs T for 
Ns/Np = 2.67, 2.13, 1.67, 1.33, 1.07, 0.94, 0.8, 0.7, 0.6, 0.5, 0.45, 
and 0.4, from top to bottom. 

from their interactions with phonons or small scattering 
sites. The equation of motion for a electron i is 



77V = f, = - vu{nj) + f/ + f/ + f,, 



(1) 



Here 77 = 1 is the damping constant and Ui{r) = —q^/r, 
with = 1, is the electron-electron interaction potential, 
treated as in [19]. The term f/ comes from the Np ran- 
domly spaced defect sites modeled as parabolic traps of 
radius Vp = 0.2 and strength fp = 1.0. The thermal noise 
arises from random Langcvin kicks with < f'^{t) >= 
and < ff{t)fj{t') >= 2rikBT6ijS{t - t'). The driving 
term = f^x comes from an applied voltage, and we 
take fd = 0.1. We start at a high temperature where 
the charges are diffusing rapidly and cool to a lower tem- 
perature. We then wait for 10^ simulation time steps 
to reduce transient effects before appling the drive and 
measuring the average velocity v of the electrons, which is 
proportional to the conduction or resistance. We do this 
for a series of electron densities at fixed disorder strength. 
We have considered samples with constant pin densities 
rip for different system sizes such that Np ranges from 
317 to 1200. 

We first consider the average electron velocity clS Qti 
function of temperature and charge density for a sys- 
tem with a fixed Np = 1200. In Fig. 1 wc show a se- 
ries of conductance curves vs T for charge density varied 
over nearly an order of magnitude, 0.4 < Ns/Np < 2.67. 
For high Ns/Np > 0.6 the conductance is finite down 
to T = 0, while for Ng/Np < 0.6, the electron velocity 
drops to zero within our resolution, indicating that all 
the electrons arc strongly pinned in an insulating phase. 
As Ng/Np increases above 1.0, the downward curvature 
of V at low T decreases. These curves appear very similar 
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FIG. 2. The relative velocity fluctuations Sv vs time for the 
system in Fig. 1 for T = 0.09 at Ns/Np = 0.5, 0.7, 1.05, 1.64, 
and 2.67, from bottom to top. The curves have been shifted 
up for clarity. 

to those typically observed in 2D MIT studies [3]. One 

difference is that we do not find a charge density n"^ 
above which the slope of the velocities turns up slightly 
at low T, as in the experiments. This may be due to 
the fact that in our model we do not directly include 
phonons. In the experimental regime of interest to us 
here, the large noise increases occur at charge densities 
Ug < n"P, where the velocity curves bend down at low 
temperature. 

We next consider the relative fluctuations in the ve- 
locities, Sv{t) = {v{t)— < V >) for varied Ng/Np at a 
fixed T — 0.09 for the system in Fig. 1. This analysis is 
similar to that performed in experiments [3-5] . In Fig. 2 
we show the time traces of the relative velocity fluctua- 
tions for Ns/Np = 0.5,0.7, 1.05, 1.64 and 2.67. Here the 
fluctuations increase as Ns drops, in agreement with the 
experiments [3]. For Ng/Np < 0.44, the system is pinned 
and there arc no fluctuations. It is possible that, over 
a longer time interval such as that accessible experimen- 
tally, there would be even larger fluctuations at these 
small Ns values; however, this is beyond the time scale 
we can access with simiilations. A similar series of time 
traces can be obtained for 5v at fixed Ns/Np for increas- 
ing temperature (not shown). Here the fluctuations are 
reduced at the higher temperatures, in agreement with 
experiments [3]. 

From the fluctuations dv we measure the power spec- 
trum 



6v{t)e- 



(2) 



In Fig. 3 we plot S{u) for two different charge densities. 
At Ng/Np = 0.5 (upper curve), the spectrum shows a 
1//" characteristic with a = 1.37 over a few orders of 
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FIG. 3. The power spectra S{i') for the velocity fluctuar 
tions Sv{t) at Ns/Np = 0.5 (upper curve) and Ng/Np = 1.67 
(lower curve). The solid line has slope a = 1.37. 

magnitude in the frequency. In contrast, for Ng/Np = 
1.67 (lower curve), the noise power at lower frequencies is 
considerably reduced and the spectra is white with a fa 0. 
We not that it is the lower frequcies which will be most 
readly accesible in experiment. For fixed Ns/Np = 0.5, 
we find that the power spectrum becomes white upon in- 
creasing T. We note that our results differ quantitatively 
from the experimental noise measurements [3] which find 
1//" noise with a = 1.0 in the metallic phase anda = 1.8 
in the glassy regime. Our exponent a = 1.37 is close to 
the Q = 1.3 found in the glassy regime in other experi- 
ments [6], where a = 0.75 in the metallic regime. It is 
possible that the exponents are not universal but depend 
on the details of the disorder strength; nevertheless, our 
results are in qualitative agreement with the experiments. 

In Fig. 4(a) we show the noise power So integrated over 
the first octave vs Ns/Np for a fixed T = 0.09. The noise 
power increases by four orders of magnitude as Ns/Np 
is reduced. At low Ng/Np, the noise power decreases 
almost exponentially with charge density and begins to 
saturate at high Ns/Np. Both these observations are in 
agreement with the experimental results [4,5]. In the in- 
set of Fig. 4(a) we plot the noise spectrum exponent a vs 
Ns/Np. A large increase in a occurs near Ng/Np = 0.7. 
A similar sharp increase in the exponent is also observed 
in experiments [4 6] as a function of charge density and 
has been interpreted as the glassy freezing transition. In 
Fig. 4(b) we plot So vs T for fixed Ns/Np = 0.5. Here 
the noise power drops exponentially over four orders of 
magnitude with increasing T, which is in agreement with 
the experiments [4,5]. We note that most single elec- 
tron models predict an increase in the noise power with 
temperature [7,8]. Other models in the hopping regime 
predict a noise power that is independent of T or has a 
power law decrease with T [20] . These discrepancies 
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FIG. 4. (a) The integrated noise power So vs Ns/Np for 
T = 0.09. Inset: The power spectrum exponent a vs N^/Np 
for T = 0.09. (b) So vs T for N^/Np = 0.5. Inset: a vs T for 
fixed Ns/Np = 0.5. 

suggest that the noise in the experiment is not due to 
single electron hopping events, but is instead caused by 
correlated electron motions. In the inset of Fig. 4(b) we 
plot a vs T, showing that a sharp increase in a occurs 
near T = 0.125 at the onset of the glassy freezing. 

We have also measured the non-Gaussian nature of 
the noise. At high Ns and T the noise fluctuations are 
Gaussian; however, in the regions of high noise power 
we find non-Gaussian noise fluctuations with a skewed 
distribution. Experiments have also found evidence for 
non-Gaussian fluctuations in the glassy regimes [5]. 

Next we show evidence that the large noise is due to 
correlated regions of string like electron flow, and that 
within these regions the electrons move in ID or quasi 
ID channels. Because of the reduced dimensionality the 
electron motion is more correlated. In Fig. 5(a) we show 
the trajectories of the electrons for a fixed period of time 
for a system with T = 0.09 at Ns/Np = 1.67. Here the 
electrons can flow freely throughout the sample, although 
there are some areas where electrons become temporarily 
trapped by a defect site. In Fig. 5(b) at Ns/Np = 1.37, 
where the noise power is larger than in the system shown 
in Fig. 5(a), larger pinned regions appear and the electron 
motion consists of a mixture of 2D and ID regions. If the 
trajectories are followed over longer times. 
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FIG. 5. Electron trajectories for a fixed period of time for 
fixed T = 0.09 at (a)Ar„/iVp = 1.67, (b) 1.37, (c)0.5, and (d) 
0.3. 



motion occurs throughout the entire sample. In Fig. 5(c) 
at Ns/Np = 0.5, where the noise power and a are both 
maximum, the electron motion occurs mostly in the form 
of ID channels that percolate through the sample. There 
are also regions where the electron motion occm-s in small 
rings. The channel structures change very slowly with 
time, with a channel occasionally shutting off while an- 
other emerges elsewhere. It is the intermittent opening 
of the ID channels which gives rise to the large noise fluc- 
tuations in this regime. When a percolating ID channel 
opens, all the electrons in that channel move in a corre- 
lated fashion leading to a large increase in the conduc- 
tion. Conversely, if a percolating channel closes all the 
electrons in that channel cease to move. It is well known 
that fluctuations in ID are much more strongly enhanced 
than in 2D. As T or Ng is increased, the motion becomes 
increasingly 2D in nature and the strong correlations of 
the electron motion are lost. We also note that the ap- 
pearance of string like motions in the large noise regions 
is consistent with studies in glassy systems, where dy- 
namical heterogeneities in the form of ID stringlike mo- 
tions of particles have been observed in conjunction with 
large noise [15]. In Fig. 5(d) at N^/Np = 0.3, deep in the 
insulating regime, there are no channels. Instead, the in- 
frequent motion of electrons occurs only by small jumps 
from defect to defect. 

In conclusion, we have presented a simple model for 
the glassy freezing of interacting electrons in 2D with 
random disorder. For high electron density or high tem- 
peratures, the electrons form a 2D liquid state and we 



find low conduction noise power with a white spectra. 
As the density of the electrons is lowered for fixed tem- 
perature, or conversely, as the temperature is lowered for 
fixed low electron density, there is a crossover to a 1//" 
noise with large low frequency power and a = 1.37. In 
this glassy regime, the electrons move in ID intermittent 
stringlike paths which percolate throughout the sample. 
Similar stringlike motions are also observed in other glass 
forming systems. For low electron density, all the elec- 
trons are frozen by the defect sites and the motion occurs 
only by single electron hopping events. We find that the 
noise power decreases exponentially with temperature, 
in agreement with experiment. Many of our results are 
in qualitative agreement with recent experiments on 2D 
electron systems near the metal insulator transition. 

This work was supported by the US DoE under Con- 
tract No. W-7405-ENG-36. 
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